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Abstract 


Analysis  of  energy  conversion  in  laser  propulsion  is  reported  and  compared  to  experimental  studies  of  a 
laboratory  scale  propulsion  device  that  absorbs  laser  energy  and  converts  that  energy  to  propellant  kinetic  energy. 
The  propellants  studied  were  air  and  Delrin,  a  solid  with  the  composition  of  formaldehyde  [H2CO]  that  vaporizes 
cleanly  upon  laser  irradiation  The  Myrabo  Laser  Lightcraft  (MLL)  was  studied.  It  incorporates  an  inverted 
parabolic  reflector  that  focuses  laser  energy  into  a  toroidal  volume  where  it  is  absorbed  by  a  unit  of  propellant  mass 
that  is  subsequently  expanded  in  the  geometry  of  the  plug  nozzle  aerospike.  The  results  showed  that  between  30  and 
50%  of  the  incident  laser  energy  is  converted  to  propellant  kinetic  energy.  This  overall  absorption/expansion 
efficiency  was  examined  in  terms  of  the  thermodynamic  predictions  of  conversion  of  propellant  internal  energy  to 
propellant  kinetic  energy.  Expansion  of  a  propellant  mass  that  was  heated  at  constant  volume  was  examined  under 
conditions  where  either  chemical  equilibrium  or  frozen  composition  was  maintained.  For  expansion  with  an 
effective  area  ratio  of  -  4,  which  is  appropriate  for  the  MLL,  a  maximum  of  25  to  50%  of  the  internal  energy  is 
predicted  to  be  convertible  to  propellant  kinetic  energy,  based  on  the  minimization  of  the  entropy  gain  of  the 
blowdown  process.  With  the  small  effective  area  ratio  —  4,  equihbrium  expansion  was  only  slightly  more  efficient 
than  frozen  expansion.  Heating  of  propellant  to  highly  ionized  states  resulted  in  lower  efficiency  energy  conversion 
but  higher  exit  velocity.  The  thermodynamic  limitations  are  illustrated  by  process  representations  of  blowdown  in 
the  Mollier  plane  (enthalpy  vs  entropy  diagram  for  air).  The  analysis  captures  the  equation  of  state  of  the  partially 
ionized  propellant  under  conditions  of  chemical  equihbrium 

Nomenclature  (in  order  of  use) 

F  Force  or  thrust,  N  =  kg  m/s2, 

m  mass  of  rocket,  kg. 

t  time,  s 

ve  exit  velocity  of  propellant  relative  to  the  rocket,  in  the  rocket  frame  of  reference,  ml s. 
dm/rit  incremental  mass  change  of  rocket,  instantaneous  propellant  mass  flow  rate,^  kg/s. 
dv/dt  incremental  velocity  change  of  rocket  in  the  inertial  frame  of  reference,  m/s'. 

Vj  initial  velocity  of  rocket  in  inertial  frame,  m/s. 
vf  final  velocity  of  rocket  in  inertial  flame,  m/s 

mi  initial  mass  of  rocket,  kg. 
mf  final  mass  of  rocket,  kg. 

mrnif  mass  of  propellant,  kg. 
f  mass  fraction  for  a  rocket  mission,  f  =  m^m;. 

x  ratio  of  vf  -Vi  to  ve. 
y  ratio  of  Vi  to  v,. 

r|i  mass  efficiency  of  mission  =  mass  of  vehicle/mass  of  propellant. 

r\2  momentum  efficiency  of  mission  =  momentum  change  of  vehicle/momentum  of  propellant 

r|3  energy  efficiency  of  mission  -  kinetic  energy  change  of  mission/kinetic  energy  of  propellant 
Ci  mass  coupling  coefficient  =  mass  of  vehicle/kinetic  energy  of  propellant 
C2  momentum  coupling  coefficient  =  momentum  change  of  vehicle/kinetic  energy  of  propellant 
C3  energy  coupling  coefficient  =  kinetic  energy  change  of  vehicle/kinetic  energy  of  propellant 
Ki  mass  cost  function  »  kinetic  energy  of  propellant/mass  of  vehicle 
K2  momentum  cost  function  =  kinetic  energy  of  propeilant/momentnm  change  of  vehicle. 

K3  energy  cost  function  =  kinetic  energy  of  propellant/kinetic  energy  change  of  vehicle. 

Sutton  Sutton's  “propulsive  efficiency”,  r[3uttxm  =  yV(  1-y2). 

M  mass  of  Laser  Lightcraft,  kg. 
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h  altitude  achieved  in  laser  powered  flight,  m. 
g  acceleration  of  gravity,  9.806  m/s2. 

El  laser  energy  incident  on  propellant  per  pulse,  J/pulse. 

a)  frequency  of  laser,  pulse/s. 

T/W  thrust  to  weight  ratio,  T/W  =  CzEi,co/Mg. 

PL  power  of  laser,  J/s,  PL  =  cdEl. 

a  efficiency  of  conversion  of  propellant  internal  energy  to  propellant  kinetic  energy. 

P  fraction  of  incident  laser  energy  absorbed  by  propellant 

I  impulse,  Ns  =  kg  m/s. 

C2(exp)  experimental  momentum  coupling  coefficient,  Ns/J  =  s/m. 
nipen  total  mass  of  pendulum  and  test  article,  kg. 

1  distance  between  pendulum  pivot  and  center  of  gravity,  m. 

0  angle  of  deflection  of  pendulum,  degrees. 

T.  Standard  impulse  imparted  to  pendulum  by  MBS  traceable  hammer  accelerometer,  Ns 

Tr-._, —  Impulse  measured  with  ballistic  pendulum.  Ns. 
me  mass  ofDelrin  ablated  in  a  single  laser  shot,  kg. 

u  specific  internal  energy  of  laser  heated  propellant,  J/kg. 

u,  specific  internal  energy  of  propellant  at  the  exit  of  the  rocket  after  isentropic  expansion. 

Introduction 

Analysis  of  the  limitations  of  laser  propulsion  and  optimization  of  the  mission  parameters  that  are  imposed 
by  Newton’s  second  law  was  published  independently  and  nearly  simultaneously  25  years  ago  by  Moeckel  (1975) 
and  Lo  (1976).  Phipps,  Reilly  and  Campbell  (2000, 2001)  recently  cited  Moeckel’s  paper  in  their  comprehensive 
analysis  of  the  Earth  to  LEO  mission  powered  by  laser  propulsion.  Here  we  show  the  elementary  relationships 
between  the  mission  parameters  and  three  mission  efficiencies,  the  mass,  momentum  and  energy  efficiency,  and 
their  associated  coupling  coefficients  and  cost  functions.  The  efficiency  of  conversion  of  laser  energy  to  propellant 
kinetic  energy,  based  on  various  ballistic  pendulum  and  flight  experiments  with  Myrabo  Laser  Lightcrafi,  MLL 
[Messitt,  Myrabo,  and  Mead  (2000);  Mead,  Squires,  Beairsto,  and  Thurston  (2000)],  is  reported. 

The  Rocket  Equation 

The  Rocket  Equation  has  been  known  for  around  one  century,  having  been  first  contemplated  by 
Tsiolkovski,  Oberth  and  Goddard.  It  has  also  been  referred  to  as  the  mass  ratio  equation  [Bussard  and  DeLauer 
(1958),  Hill  and  Petersen  (1965)]  and  the  basic  Ziolkowsky  equation  [Lo,  (1976)].  The  thrust,  F,  that  results  from 
expulsion  of  matter  from  a. vehicle  is  expressed  by  Newton’s  second  law,  as 

(1)  F  =  d(mv,)/dt, 

where  mv,  is  the  momentum  of  the  jet  exhaust  in  the  vehicle  frame  of  reference,  [Corliss,  I960],  For  the  case  where 

v,  is  constant,  differentiation  of  Equation  (1)  produces 

(2)  F  =  v,  dm/dL 

Figure  1  shows  the  rocket  parameter  space  based  on  Equation  (2)  spanning  14  orders  of  magnitude  in  thrust  The 
power  of  the  rocket  jet  is  P  =  Vi  dm/dt  v.2  =  i4Fve.  Specification  of  any  two  of  the  four  rocket  parameters  of  interest 
(thrust  exit  velocity,  power,  and  propellant  mass  flow  rate),  fixes  the  other  two  parameters.  Equation  (1)  may  also 
be  used  to  an  average  exit  velocity  for  rockets  where  v,  is  not  constant  such  as  blowdown  from  a  fixed 
volume.  Thus,  v^„  is  the  total  impulse  (the  integral  of  Fdt  over  the  blowdown  time)  divided  by  the  total  mass  (the 
integral  of  dm  over  the  blowdown  time),  so  that  v^,  is  the  mass  weighted  average  exit  velocity. 

The  Rocket  Equation  results  from  a  balance  of  the  force  exerted  by  the  propellant  on  the  vehicle  and  the 
motion  of  the  vehicle  under  the  influence  of  the  propulsive  force  as  required  by  Newton’s  second  law.  Thus,  in  the 
absence  of  other  forces,  such  as  body  (gravitational)  force  and  drag  force. 
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(3 )  F  =  v,dm/dt  =  -  mdv/dL 


where  v  is  the  vehicle  velocity  in  the  inertial  frame  of  reference,  Le.,  the  velocity  relative  to  a  fixed  point  is  space. 
Elimination  of  time!  in  Equation  (1)  yields  the  expression  for  conservation  of  momentum,  mdv  =  -v.dm,  which  may 
be  integrated  between  the  limits  of  initial  and  final  mission  velocity  (v*  and  vf)  and  mass  (mi  and  m*)  to  produce  the 
Rocket  Equation, 

(4)  f  =  m/mi  =  exp-{vrVi)/v,  =  exp(-x). 

Mission  Efficiencies 

The  mass,  momentum  and  energy  efficiencies  for  the  mission  defined  by  Vf,  vt  and  v0  follow  naturally  from 
the  Rocket  Equation,  viz., 

(5)  rii  =  nif/Cmi-mf)  =  f/(l-f)  =  l/(ex-l), 

(6)  t|2  =  (mfVrniiVi)/(nii-mt)ve  =  fr(l-f)vA.  -  l/(l-f)Vi/v.  =  x/(ex-l)  -  y, 

(7)  r|3  =  (mfV^-mjVi^/Cmrmdv.^l-flCVf/v,)2  -  l/(l-f)(vj/ve)2  =  (x2+2xy)/(ex-l)  -  y2 , 
where  y  =  v/v,. 


For  missions  with  zero  initial  velocity,  t|i  -  l/(ex-l),  t|2  —  x/(ex-l),  and  ri3  —  x2/(ex- 1).  Figures  2a,  2b,  and 
2c  show  three  representations  of  the  change  in  the  propulsion  efficiencies  as  the  mission  progresses  from  f  =  1 
toward  f=0.  Inthef=  1  limit,  m  is  infinite,  r|2is  unity,  and  q3  is  zero.  Also,  in  this  limit,  r)3  =  1/pi  =  x.  Whenf= 
1/e  =  0.367879,  pi  =  t|2  =  r|3  =  0.581976.  A  maximum  in  the  energy  efficiency,  r|3  =  0.647635,  occurs  atf = 

0.203 179,  Vf/ve  =  1.59367.  Since  these  limits  are  fundamental  mathematical  constants  derived  from  e  =  2.71828, 
they  may  be  written  with  as  many  significant,  figures  as  desired.  The  limits  may  be  established  by  expanding  the 
exponential  in  a  Maclaurin  series,  viz.,  f1  =  ex  =  1  +  x  +  xV2!  +  ...,  and  truncating  after  the  second  term  since  x  « 

1.  The  earliest  publication  of  the  maximum  in  t|3,  for  the  case  of  zero  initial  velocity,  appears  to  be  in  Lo’s  1976 
paper,  “Propulsion  by  laser  energy  transmission.”  Phipps,  et  aL  (2000,  2001),  optimized  their  Earth  to  LEO  mission 
with  Vf/v„  =  1.59367.  The  Figures  show  that  the  momentum  efficiency,  initially  unity,  falls  below  0.95  only  after  the 
mass  fraction  decreases  below  about  0.9.  Missions  with  f  >  0.99,  which  includes  short  Laser  Lightcraft  flights  and 
bench  top  measurements  of  momentum  (Le.,  impulse,  with  single  laser  shots),  have  t|2  >  0.995. 

Upper  stage  missions  have  initial  velocities  greater  than  zero.  The  mass  efficiency  is  independent  of  the 
initial  velocity,  but,  as  shown  in  Figure  3a,  momentum  efficiencies  decrease  as  V|  increases,  and,  as  shown  in  Figure 
3b,  energy  efficiencies  increase  as  Vi  increases.  Figure  3b  shows  that  as  Vj/v,  increases  from  zero  to  unity,  the 
maxima  in  t)3  increase  toward  unity  and  shift  toward  larger  f  [or  toward  smaller  (vi-vi)/v,].  Thus,  the  missions  of 
upper  stage  rockets  may  be  accomplished  with  higher  energy  efficiency.  For  example,  if  a  second  stage  has  a 
propellant  with  exit  velocity  v^  equal  to  its  initial  velocity  Vq,  ti3  is  unity  when  f  =  1,  at  the  instant  of  firing  the 
upper  stage.  As  the  mission  progresses,  v  increases,  f  falls  below  unity,  and  q3  fells  and  eventually  becomes 
negative  when  f  <  0.0813.  At  f  =  0.0813  the  kinetic  energy  of  the  vehicle  is  actually  unchanged  from  what  it  was 
at  the  beginning  of  the  mission  but  the  vehicle  velocity  is  substantially  increased  to  a  value  around  2.5 1  times  the 
second  stage  propellant  exit  velocity,  v^. 

Nearly  forty  years  ago  Sutton  (1966)  introduced  the  notion  of  “propulsive  efficiency”,  which  we  denote 
Ttsunonj  to  characterize  the  affect  of  Vi  on  the  initial  value  of  ri3,  i.e.,  =  limit  as  f  approaches  umty  of  p3:  T|sunm  = 

y /( 1  -y2)  where  y  =  Vi/v,.  Sutton  described  perfect  kinetic  energy  transfer  from  the  jet  to  the  vehicle  as  the  condition 
achieved  when  Vi  =  v,  so  that  propellant  is  deposited  with  zero  velocity  in  the  inertial  frame  of  reference. 
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Coupling  Coefficients  and  Cost  Functions 


Mass,  momentum  and  energy  coupling  coefficients  (Q,  C2  and  C3, 
coupling  of  vehicle  parameters  to  the  kinetic  energy  of  the  propellant  They 
energy  of  the  vehicle  per  unit  of  propellant  kinetic  energy,  and  are  simply  related  to  the  propulsion  efficiencies. 


respectively)  indicate  the  degree  of 
are  defined  as  the  mass,  momentum  and 


(8)  Ci  =  t|i/v,2. 


(9)  C2  =  2  tia/v. 


(10)  C3  =  T|3 

Figure  4  shows  the  dependence  of  the  momentum  coupling  coefficient,  C2,  on  v.,  vf  and  f  for  missions 
where  the  initial  velocity  v,  =  0.  Figure  4  may  be  used  to  show  the  connection  between  experimental  measurements 

of  C2  and  low  altitude  flights,  where  f  is  near  unity,  and  Earth  to  LEO  missions,  where  f  for  maximum  energy 

propulsion  efficiency  is  0.203.  Figure  4  also  shows  that  in  the  f  =  1  limit,  C2  =  2/v„  which  is  the  condition  achiev 
in  the  experimental  work  reported  here. 

The  energy  coupling  coefficient,  C3,  and  the  energy  efficiency,  ri3,  are  identical  and  depend  only  on  the 
mass  fraction  and  initial  velocity  as  has  been  previously  shown  in  Figure  3b. 


The  mission  cost  functions  express  mission  costs  in  terms  of  investment  of  propellant  kinetic  energy  ger-® 
unit  of  vehicle  mass,  momentum,  and  energy  for  a  mission  specified  by  vf,  Vi,  and  ve.  The  cost  functions  are 
reciprocals  of  their  respective  coupling  coefficients: 


(11)  Ki  =  1/Ci, 

(12)  K2  =  1/C2, 


(13)  K3  =  1/C3. 

Phipps,  etaL  (2000,  2001)  used  the  Ki  cost  function  (kinetic  energy  of  propellant  per  unit  of  propelled 
mass)  in  their  Earth  to  LEO  mission  analysis  to  illustrate  the  optimization  of  mission  parameters.  The  first  use  01 
appears  to  be  in  Moeckel’s  1976  paper,  “Optimum  exhaust  velocity  for  laser  driven  rockets.’  Figure  5  shows  the 
parameter  relationships  of  Ki  for  missions  with  zero  initial  velocity.  The  energy  cost  is  defined  as  propellantionetic 
energy  investment  per  unit  of  propelled  mass.  The  Figure  shows  that  each  mission,  specified  by  vf,  has  an  optimum 
Ve,  such  that  Ki  is  minimized  for  v,/v.  =  1.5937.  It  should  be  noted  however  that  the  minima  m  K*  for  the  lines  of 
constant  v,  are  quite  shallow,  and  that  expenditure  of  only  slightly  more  energy  may  be  desirable  because  ot  me 
attendant  relatively  larger  gain  in  mass  traction. 


The  analysis  of  a  laser  powered  Earth  to  LEO  mission  by  Phipps,  et  al  (2000, 2001)  included  gravity  and 
drag  losses  and  optimization  of  mission  time  to  minimize  these  losses.  For  their  vehicle,  they  concluded  mat  the 
required  effective  vt  was  around  10,000  m/s,  so  that  gravity  and  drag  losses  amounted  to  an  equivalent  velocity  loss 
of  about  2,000  m/s  over  a  flight  time  ~  800  s.  As  shown  in  Figure  4,  a  propellant  with  exit  velocity  of  about  6000 
m/s  and  an  overall  mission  momentum  coupling  coefficient  of  about  1.2  x  10  Ns/J  would  accomplish  the  single  ^ 
stage  mission  to  LEO  with  a  mass  fraction  of  0.203,  which  translates  to  propulsion  efficiencies  of  m  -  0.255,  q2  - 
0.406,  and  Ti3  =  0.648. 


Flights  of  Laser  Liehtcraft 

Knowledge  of  the  momentum  coupling  coefficient  alone,  Cj,  allows  prediction  of  the  altitude  as  a  function 
of  time  in  short  flight  experiments.  For  constant  C2  during  the  flight,  and  neglecting  drag,  the  balance  of  forces  on 
the  vehicle  of  mass  M  may  be  written 

(14)  F  =  M(d2h/dt2-g)  =  C2ELffl, 
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where  h  is  altitude,  to  is  laser  frequency,  and  El  is  the  energy  of  one  laser  pulse  incident  on  the  propellant.  Equation 

(14)  is  appropriate  for  air,  where  M  remains  constant,  and  for  Delrin  propellant  in  flight  experiments  that  utilized 
very  small  propellant  mass  so  that  M  is  effectively  constant,  f  >  0.99,  and  r\2  ~  1.  Integration  of  Equation  (14)  from 
time  zero,  where  h  =  0,  dh/dt  =  0,  d2h/dr  =  0,  yields  altitude  as  a  function  of  time,  coupling  coefficient,  C2,  laser 
energy  incident  on  propellant  per  pulse,  El,  and  laser  pulse  frequency,  a : 

(15)  h  =  t2/2(C2ELffl/M-g). 

The  thrust  to  weight  ratio  alone,  T/W,  determines  the  flight  characteristics,  T/W  =  C2ELco/Mg.  Since  TAV  >  1,  is 
required  for  flight,  C2  >  Mg/ELco  is  necessary,  where  ELco  is  the  effective  laser  power,  PL.  Figure  6  shows  predicted 
flight  paths  for  several  values  of  C2  that  span  measured  values  in  single-shot  bench  top  experiments  with  air  and 
Delrin.  The  flight  paths  are  shown  for  a  0.04  kg  vehicle  propelled  by  10  kW  of  laser  power  at  to  =  25  pulses  per 
second,  EL  =  400  J/pulse.  Figure  7  shows  a  time  exposure  photograph  of  a  MLL  flight  A  new  altitude  record,  h  ~ 
80  meters  was  recently  achieved  with  air  propellant  (Myrabo,  October  2000).  Obviously,  by  comparing  predicted 
flight  paths  to  time  exposure  photography  or  video  tape  recordings,  coupling  coefficients  and  their  variation  during 
the  flight  may  be  determined. 

Energy  Conversion  Efficiency  of  Laser  Propulsion 

The  efficiency  of  conversion  of  laser  energy  to  propellant  kinetic  energy  may  be  defined  by  energy 
conservation, 

(16)  propellant  kinetic  energy  =  Vifmi-mdv,,2  =  ap  EL, 

where  El  is  the  laser  energy  incident  on  the  propellant,  p  is  the  fraction  of  energy  absorbed  into  a  mass  m-m.-  of 
propellant,  and  a  is  the  fraction  of  the  absorbed  energy  that  is  converted  into  propellant  kinetic  energy  by  expansion. 
The  value  of  ap  may  be  determined  by  experimental  measurement  of  the  impulse  imparted  to  a  test  article  when  a 
known  amount  of  propellant  mass  is  ejected  by  a  laser  pulse  of  energy  El.  Since 

(17)  v.  =  I/(mrmf), 

(18)  ap  =  I2/2EL(mi-mf), 

which  may  also  be  written  in  terms  of  the  commonly  reported  experimental  coupling  coefficient,  C2(exp)  =  I/El, 

(19)  ap  =  C2(exp)  I/2(mrmf). 

Thus,  to  determine  ap,  the  propellant  mass  associated  with  the  measured  impulse  and  laser  pulse  energy 
must  also  be  measured. 

The  kinetic  energy  transferred  to  a  vehicle  (in  free  flight  or  mounted  on  a  ballistic  pendulum)  is  equivalent 
to  the  potential  energy  it  achieves,  mfgh.  In  experiments  with  f  >  0.999,  this  energy  is  a  very  small  fraction  of  the 
total  jet  kinetic  energy,  as  shown  by  Equation  (7)  and  Figure  2.  However,  if  the  potential  energy  and  propellant 
mass  are  measured,  a  second  independent  determination  of  ap  is  possible: 

(20)  aP  =  (nVmi-nif)  ntigh/EL. 

Thus,  two  methods  for  measurement  of  ap  are  available,  both  of  which  require  measurement  of  propellant  mass,  mr 
nif.  In  the  impulse  measurement,  error  in  the  ap  determination  is  magnified  because  ap  is  proportional  to  I2, 
whereas  in  the  potential  energy  measurement,  the  ap  error  is  proportional  to  the  measured  h  of  a  flight  experiment 
or  the  measured  angle  of  deflection  in  a  ballistic  pendulum  experiment  (vide  infra). 
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Experimental 

A  ballistic  pendulum  was  used  to  measure  the  impulse  imparted  to  test  articles  by  absorption  of  a  pulse  of 
laser  energy  followed  by  expansion  of  the  heated  propellant.  The  10  kW  pulsed  CO2  laser  maintained  by  the  High 
Energy  Laser  System  Test  Facility  (HELSTF)  at  White  Sands  Missile  Range,  NM  was  used.  The  laser  produces 
pulse  energies  of  10.6  pm  light  up  to-  500  J  per  pulse  over  10  to  30  ps  at  repetition  rates  up  to  about  30  Hz.  Laser 
pulse  energy,  El,  was  measured  with  a  calorimeter  and  is  accurate  to  better  than  10%. 

The  three  test  articles  studied  were  Myrabo  Laser  Lightcraft  (MLL),  which  is  an  inverted  paraboloid 
reflector  that  also  acts  as  a  plug  nozzle,  as  illustrated  in  Figure  8.  [see  for  example,  Mead,  et  al,  (2000),  Figure  4,  for 
a  cross-sectional  view].  The  focusing  optic  was  identical  in  each  test  article.  It  focused  a  square  beam,  -  10  cm 
with  a  4  cm2  antral  square  void,  to  a  circle  of  laser  light  that  fell  on  the  surface  of  the  aluminum  or  SiC  shroud 
With  air  propellant,  the  energy  is  absorbed  into  a  volume  of  air  no  larger  than  -  25  cm3  which,  at  standard 
temperature  and  pressure,  amounts  to  -  25  mg  of  air  maximum  This  upper  limit  to  the  mass  of  air  propellant  is 
based  on  the  geometry  of  the  MLL  model  200-3/4  geometry.  With  Delrin  propellant,  the  mass  of  ablated  propellant 
was  determined  by  weighing  before  and  after  each  series  of  ablation  experiments  at  fixed  EL.  Table  I  summarizes 
the  test  article  masses,  pendulum  total  masses,  m^  and  distances,  1,  between  the  pivot  and  center  of  gravity. 


Table  I.  Experimental  test  articles  and  pendulum  parameters,  MLL  model  200-3/4. 


article 

article  mass, 
kg 

pendulum  mass, 
HU  kg 

cofg,l, 

m 

propellant 

bare  pendulum 

0.7035113 

0.7035113 

0.3556 

1.  heavyweight  model, 
aluminum  shroud 

0.3140756 

1.0175869 

0.3826 

air 

2.  SiC  shroud  flight 
model 

0.1141644 

0.8176757 

0.3683 

air 

3.  SAR  (solid  ablative 
rocket)  flight  model 

0.1241798 

0.8276911 

0.3690 

Delrin  plus 
air 

The  experimental  impulse  was  obtained  from  the  measured  maximum  angle  of  deflection,  9,  of  the 
pendulum  during  the  first  half  cycle  after  imparting  the  impulse: 

(21)  I  =  mfM1[gl(l-cos9)]1/2  -  m^xO  [gl]1/2/720, 

where  m^  is  the  pendulum  plus  test  article  mass,  g  is  the  gravitational  constant,  9.806  m  s'2,  and  1  is  the  distance 
between  pendulum  pivot  and  center  of  gravity  [Thomas  (1953)].  The  approximation  to  Equation  (21),  which  derives 
from  the  Maclaurin  expansion  for  cos9,  is  accurate  to  better  than  one  part  in  ten  thousand  when  9  <  10  degrees. 
Equation  (21)  derives  from  the  equivalence  between  the  initial  kinetic  energy  and  the  potential  energy  when  the 
pendulum  is  at  maximum  deflection,  !4  mpenvpal2==mpmgl(l-cos9)  under  conditions  where  there  are  no  losses.  Based 
on  the  measured  damping,  the  energy  loss  was  about  0.00015  J  (0.6  %)  per  half-cycle  of  the  total  of  0.0230  J  of 
potential  energy  when  the  largest  value  of  9,  -  7  degrees,  was  measured  so  no  correction  for  this  loss  was  made  to 
the  data 

The  mass  of  propellant  was  typically  a  few  tens  of  milligrams,  which  amounts  to  amass  fraction  off- 
0.9999  to  0.99999.  Thug  whereas  momentum  from  the  jet  is  transferred  to  the  pendulum  with  unit  efficiency,  r|2  = 
1,  the  energy  of  the  jet  is  transferred  to  the  pendulum  with  efficiency  ri3  — 10"4  to  10 5  (see  Figure  2). 

For  each  of  the  test  articles,  the  pendulum  was  calibrated  with  an  NBS  traceable  impulse  hammer.  The 
hammer  imparted  impulses  ranging  from  0.005  to  0.04  N-s  in  about  15  ms.  Figure  9  compares  the  impulse  from  the 
hammer  calibrations  to  those  obtained  from  the  pendulum  The  figure  shows  that  the  ratio,  Ipendsiam/Itwmmer,  varied 
within  about  2  %  of  unity,  and  indicates  that  Ip,„duian.  is  systematically  higher.  This  is  probably  due  to  a  small 
energy  loss  in  the  hammer  calibration  that  occurs  during  the  imperfect  impact  of  the  hammer  against  the  pendulum 
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Results 


Figure  10  shows  the  mass  of  Dekin  ablated,  mo,  as  a  function  of  El  for  the  SAR  test  article,  and  Figure  1 1 
shows  the  exit  velocity  based  on  the  measured  impulse,  v, =  I/uid-  Since  a  small  amount  of  air  could  also  become 
involved  in  the  propellant  expansion,  the  v,  values  represented  here  are  upper  limits,  mo  was  measured  to  an 
accuracy  of  a  few  tenths  of  milligrams  by  difference  weighing  with  a  Mettler  balance. 

Figure  12  shows  the  dependence  of  the  overall  energy  conversion  efficiency,  a(3,  onEL  for  the  SAR  test 
article.  Since  these  values  of  a.p  are  based  on  the  measured  impulse,  they  do  not  account  for  possibly  significant 
pendulum  losses  that  may  occur  when  the  laser  thrust  profile  [a  few  hundred  microseconds  long,  (Jones,  2001)]  is 
imposed  on  the  Lighlcraft  These  losses  are  nearly  absent  in  the  hammer  calibrations  where  the  thrust  profile  is 
stretched  over  a  one  hundred  fold  longer  time  scale.  The  high  instantaneous  thrust  of  the  laser  initiated  detonation 
may  lead  to  significant  loss  in  the  somewhat  loose  moving  parts  of  our  ballistic  pendulum.  If  these  losses  were 
quantified  oc(3  would  increase.  It  should  be  pointed  out  that  this  series  of  pendulum  measurements  was  immediately 
followed  by  irtentiral  experiments  where  impulse  was  measured  with  a  fast  piezo  electric  transducer  capable  of 
better  than  1  ms  time  resolution.  The  piezo  electric  measurements  (Jones,  2001)  showed  that  the  integrated  impulse 
were  nearly  two  times  the  impulse  measured  with  the  pendulum  This  means,  since  ap  goes  as  the  square  of  the 
impulse,  that  ap  may  be  significantly  larger  than  has  been  derived  here  from  pendulum  measurements  without 
consideration  of  losses. 

Figure  13  summarizes  selected  momentum  coupling  coefficients,  C2,  obtained  for  air  and  Delrin 
propellants  over  the  course  of  several  experimental  test  periods  at  WSMR  during  the  past  few  years  with  our 
ballistic  pendulum  There  appears  to  be  a  strong  dependence  of  C2  on  the  quality  of  the  beam  as  between  a  tightly 
focused  beam  that  produces  lower  C2  with  air  propellant  than  a  loosely  focused  beam  This  may  be  due  to  the  tight 
beam  heating  a  gmaiiw  mass  of  air  to  a  higher  energy  than  the  more  diffuse  loosely  focused  beam.  Although  the 
exit  velocity  is  higher  in  the  tight  beam  case,  the  total  impulse  is  lower  because  the  heated  mass  is  lower. 
Comparison  of  air  to  Delrin  propellant  with  a  tightly  focused  beam  shows  that  with  Delrin  the  measured  C2  values 
are  two  to  three  times  larger,  probably  because  a  larger  mass  of  Delrin  is  ablated  than  is  heated  in  the  case  of  air. 

Discussion  of  Thermodynamic  Considerations 

Figure  14  shows  the  chemical  equilibrium  Mollier  diagram  for  air  up  to  24,000  K.  Figure  14  is  based  on 
the  database  maintained  at  NASA/Glenn  (McBride  and  Gordon  (1996)  which  is  certified  accurate  up  to  20,000  K 
and  based  on  extended  9-parameter  fits  to  enthalpy,  heat  capacity,  and  entropy  of  neutral  species  and  singly  charged 
inns  Above  20,000  K  doubly  charged  ions  begin  to  contribute  but  these  are  not  included  in  the  database. 

In  Figures  15  and  16  we  show  an  abbreviated  pair  of  Mollier  diagrams  for  air  which  we  use  to  illustrate 
process  representations  of  conversion  of  internal  energy  to  kinetic  energy  and  to  compare  this  energy  conversion 
when  isentropic  expansions  occur  under  conditions  of  equilibrium  (Figure  15)  or  frozen  flow  (Figure  16).  This  is  a 
highly  simplified,  single  parameter  analysis  based  on  the  notion  of  an  absorption  volume  into  which  a  fraction,  p,  of 
the  mriHpnt  laser  energy  is  deposited.  Thus,  we  may  consider  the  range  of  internal  energy  appropriate  for  the  model 
200-3/4  MLL  that  would  be  applicable  to  laser  heated  air  in  our  experiments.  The  Figures  show  a  scale  of  specific 
internal  energy  with  p  —  1  based  on  the  three  conceivable  volumes  of  laser  heated  air  at  STP  density  (1.18  kg/m  ), 
5.4, 21.7,  and  43.4  cm3.  These  are  the  volumes  of  toroids  with  0.5, 1.0,  and  1.5  cm  diameters  that  fit  snugly  into  the 
shroud  of  the  model  200-3/4  MLL.  The  scales  show  that  a  specific  internal  energy  of  4e7  J/kg  would  be  reached 
when  about  260  J  of  laser  energy  is  deposited  into  5.4  cm3  or  when  about  1000  J  of  laser  energy  is  deposited  into 
21.7  cm3  of  STP  air.  Additionally,  the  state  with  specific  internal  energy  of  2e7  J/kg  would  be  reached  with  about 
120  J  deposited  in  5.4  cm3  or  about  500  J  in  21.7  cm3  of  STP  air.  The  Figures  also  show  three  isentropes  originating 
from  the  constant  u  =  4e7  J/kg  line  and  the  constant  u  =  2e7  J/kg  line.  These  isentropes  are  labeled  wih  circles  that 
inHiratfi  the  area  ratios  for  the  isentropic  expansion. 

Figure  15  shows  that  equilibrium  blowdown  from  the  4e7  J/kg,  1.18  kg/m3  state  with  an  area  ratio  of  4 
produces  an  exit  pressure  of  -  1  atm  and  exit  internal  energy,  u,  ~  2.8e7.  The  exit  velocity  is  thus  v„  =  [2(4e7- 
2.8e7)J/kg]lfl  =  4900  m/s  and  the  fraction  of  energy  converted  to  kinetic  energy  is  a  =  1.2e7/4e7  =  0.3.  In 
comparison,  blowdown  from  the  2e7  J/kg,  1. 18  kg/m3  state  with  an  area  ratio  of  4  also  produces  an  exit  pressure  of- 
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1  atm  and  the  exit  internal  energy,  u*  ~  1.2e7.  The  exit  velocity  is  thus  ve  =  [2(2e7-1.2e7)  J/kg] 1/2  =  4000  m/s  and  the 
fraction  of  energy  converted  to  kinetic  energy  is  ct  =  0.8e7/2e7  =  0.4.  Thus,  the  indications  are  that  laser  heating  to 
higher  internal  energy  results  in  lower  energy  conversion  efficiency  but  higher  exit  velocity. 

Comparison  of  frozen  blowdown  (Figure  16)  shows  that  frozen  flow  loss  occurs.  With  area  ratio  of  4,  from 
the  u  =  4e7  J/kg,  p  =  1.18  kg/m3  initial  state,  ve  =  4500  m/s  and  a  =  0.22  (loss  of  ~  400  m/s  in  vt  and  0.08  in  a),  and 
from  the  u  =  2e7  J/kg  state,  v.  =  3200  m/s  and  a  =  0.25  (loss  of  800  m/s  in  v,  and  0. 15  in  a).  Comparison  of  Figures 
15  and  16  shows  that  frozen  losses  increase  drastically  when  the  effective  area  ratio  increases  beyond  4. 

Figure  17  shows  the  dependence  of  a  on  area  ratio  and  initial  density  for  isentropic  blowdown  from  u  = 

4e7  J/kg  to  an  exit  pressure  determined  by  the  area  ratio.  In  Figure  17  only  the  energy  conversion  from  momentum 
thrust  is  captured  At  area  ratio  of  4  the  exit  pressure  is  -  1  atm.  With  an  initial  density  of  STP  air  (1. 18  kg/m3)  the 
Figure  shows  a  =  0.3 1.  In  Figure  18  the  efficiency  of  energy  conversion  from  both  momentum  and  pressure  thrust 
by  blowdown  to  vacuum  is  shown.  Similarly,  for  u  =  4e7  J/kg,  area  ratio  of  4,  initial  density  of  STP  air,  the  Figure 
shows  a  =  0.40,  so  that  operation  in  vacuum  produces  a  substantial  (~  25%)  increase  in  the  conversion  of  internal 
energy  to  propellant  kinetic  energy. 

Figure  19  shows  exit  velocity  from  equilibrium  expansion  of  air  heated  to  u  =  4e7  J/kg  and  initial  densities 
ranging  from  le-5  to  le2  kg/m3  as  a  fimction  of  area  ratio.  Curves  with  filled  triangles  denote  exit  velocity  that 
derives  from  momentum  thrust  alone  and  are  transformations  of  isentropes  from  the  Mollier  plane  that  originate 
from  the  constant  u  —  4e7  J/kg  line.  Curves  defined  by  squares  show  exit  velocity  for  expansion  into  vacuum  and 
include  impulse  from  momentum  and  pressure  thrust  Curves  defined  by  circles  show  exit  velocity  for  expansion  to 
1  atm  external  pressure.  These  are  truncated  at  the  point  where  the  impulse  from  pressure  thrust  becomes  negative 
due  to  over-expansion  of  the  propellant  for  the  given  area  ratio.  The  actual  average  exit  velocity  from  the  blowdown 
beginning  at  a  given  density  is  the  density  weighted  average  exit  velocity.  For  example,  the  heavy  solid  line 
represents  the  vacuum  exit  velocity  for  blowdown  from  an  initial  density  of  1  kg/m3  and  an  area  ratio  of  4.  Initially 
the  exit  velocity  is  ~  5800  m/s  and  the  nozzle  exit  density  is  ~  0.09  kg/m3.  After  90%  of  the  propellant  has  blown 
down,  the  exit  velocity  is  -  5500  m/s  and  the  nozzle  exit  density  is  -  0.009  kg/m3.  After  99%  blowdown,  the  exit 
velocity  is  —  5300  m/s  and  the  nozzle  exit  density  is  ~  0.0009  kg/m3.  At  99.9%  blowdown  exit  velocity  is  5100  m/s 
and  the  nozzle  exit  density  is  0.0.00009  kg/m3.  Thus  the  average  vacuum  exit  velocity  is  ~  5763  m/s,  is  only  -  1% 
less  than  the  initial  exit  velocity.  In  this  example,  if  expansion  to  external  pressure  of  1  atm  pressure  at  area  ratio  of 
4  is  considered,  the  initial  exit  velocity  is  5500  m/s  and  after  90%  blowdown  it  is  ~  5000  m/s.  Thus,  there  is  a 
significant  negative  pressure  thrust  involved  when  blowdown  to  1  atm  external  pressure  occurs. 

Conclusions 

Experimental  studies  of  the  200-3/4  model  Myrabo  Laser  Lightcraft  with  air  and  Delrin  propellants  heated 
by  10.6  m  radiation  from  a  C02  laser  showed  that  energy  conversion  efficiencies  of  laser  energy  to  propellant 
kinetic  energy  were  at  least  30%  and  may  be  as  high  as  50%.  This  was  found  to  be  consistent  with  highly  simplified 
analysis  of  the  thermodynamics  of  equilibrium  from  initial  states  specified  by  a  single  parameter,  the  volume  into 
which  the  laser  energy  is  absorbed.  The  measured  exit  velocity  based  on  measured  Delrin  mass  and  estimated  air 
mass  are  in  the  neighborhood  of  ~  2000  —  3000  m/s,  which  is  somewhat  less  than  predicted  by  the  thermodynamic 
model.  In  the  case  of  Delrin,  it  was  noted  that  laser  pulse  energies  in  excess  of  -  350  J/pulse  frequently  caused 
fragmentation  and  ejection  of  large  macro  chunks  of  propellant  —  500  mg,  when  clean  ablation  produced  only  about 
30  mg  of  propellant  The  relatively  low  exit  velocity  measured  in  the  case  of  Delrin,  -  2000  m/s,  is  quite  probably 
due  to  the  ablation  of  a  significant  fraction  material  as  molecules  or  assemblies  of  molecules  in  addition  to  hydrogen 
and  carbon  monoxide.  It  should  be  noted  that  beam  quality  plays  an  important  role  in  the  performance  of  the  model 
200-3/4  MLL  which  is  counterintuitive  inasmuch  as  lower  beam  quality  (energy  spread  over  a  larger  area)  produces 
higher  Cz  coupling  coefficients  in  the  bench  top  and  flight  experiments  of  the  past  Computational  fluid  dynamics 
modeling  of  the  absorption  (and  reflection)  of  laser  energy  and  expansion  of  the  formed  plasma  [Wang,  Chen,  Liu, 
Myrabo,  and  Mead  (2001)]  have  recently  been  carried  out  The  simple  analysis  presented  here  may  only  be  useful 
in  providing  upper  limitations  to  the  conversion  of  laser  energy  to  propellant  kinetic  energy  and  to  provide  a 
simplified  visualization  and  description  of  the  processes  occurring  in  blowdown  of  laser  heated  propellants  in  MLL 
devices. 
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Figure  3a.  Momentum  efficiency  for  missions  with  various  initial  velocity ,  v/ve  =  0  to  L 


Figure  3b.  Energy  efficiencies  for  missions  with  various  initial  velocity ,  v/ve  ~  0  to  1. 


2 


Mass  cost  function,  Ki(  (J  of  propellant  kinetic  energy /kg  of  propelled  mass) 


2 


1.0  0.9  0.8  0.7  0.6  0.5  0.4  0.3  0.2  0.1  0.0 

f  =  exp-(v/v8)  _ 


Figure  4.  Momentum  coupling  coefficient  for  missions  with  zero  initial  velocity  as  a  function  of  ve,  vj,  andf 
Lines  of  constant  vf  (dotted  lines)  and  ve  (solid  lines)  are  shown  at  increments  of  1000  m/s  above  2000  ms.  In  the 
limit  of  f  -  I,  C2  =  2/ve.  The  Earth  to  LEO  mission  requires  10, 000  m/s  including  gravity  and  drag  losses, 

which  is  optimized  at  f  ~~  0.2  with  vP  =  6000  m/s  and  C?  ~  0.000 1  Ns/J  f Phipps .  Reilly  and  Campbell  (2000)1. 


Figure  5.  Mass  cost  function,  Kh  as  used  by  Phipps,  Reilly  and  Campbell  (2000)  in  Earth  to  LEO  mission 
analysis.  The  value  of  K/  is  optimized for  all  missions  with  zero  initial  velocity  when /=  0.203. 
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Figure  6.  Laser  Lighter  aft  flights  with  air  with  various  coupling  coefficients  for  a  Model  200-3/4  laser 
Lightcraft,  with  mass  of  40  grams .  These  flight  patterns  are  also  valid  for  flights  with  onboard  propellant 
provided  that  the  mass  of  ablated  propellant  is  much  less  than  the  test  article  mass,  f>  0.99 .  T/W  =  1  when 
C2  -  39.2  Ns/MJ.  _ 


Figure  7.  Time  exposure  of  nighttime  flight  of 
Myrabo  Laser  Lightcraft.  The  time  between 
flashes  of  the  air  plasma  is  0.04  s. 


Figure  8.  The  Myrabo  Laser  Lightcraft  showing 
air  plasma.  The  Model  200-3/4  is  ~  0.1  m 
diameter  at  largest  circumference.  The  aluminum 
model  weighs  ~  30  g  without  Delrin.  About  10  g  of 
Delrin  was  used  in  the  Solid  Ablative  Rocket  (SAR) 
of  which  M).3  g  was  ablated  during  a  typical  flight 
with  about  100  shots .  _ 
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Figure  9 .  Comparison  of  Impulse  measured  by  ballistic  pendulum  to  impulse  measured  with  standard  hammer 
accelerometer. 


20  Oct  2000  200-3/4  SAR  with  del rin. 

124.1978  gram  200-3/4  SAR...calib  runs  31-45 

curve  fit  with  2nd  order  regression: 

mass  =  -  2.6059e-2  +  0.161*EL-  1.3297e-4*EL2 


Figure  ID.  Mass  of  Delrin  evaporated  per  laser  pulse  of  given  energy.  The  points  labeled  with  a  \ J *  were 
accompanied  by  impulse  measurements  with  a  piezoelectric  transducer  which  indicated  that  the  duration  of  the  thrust 
was  applied  over  a  time  of  several  hundred  microseconds  (Jones,  2001).  The  width  of  the  laser  pulse  was  ~  18  us. 
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Figure  1 L  Exit  velocity  of  Delrin  based  on  measured  impulse  and  measured  mass  of  ablated  Delrin. 


Figure  12 .  Fractional  conversion  of  laser  energy  to  propellant  kinetic  energy  (aj5)  with  Delrin  propellant. 
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Figure  13.  Comparison  of  coupling  coefficient  measured  for  air  and  Delrin  with  a  tightly  focused  laser  beam  and 
comparison  of  tightly  focused  and  loosely  focused  laser  beam  with  air  propellant.  The  lower  quality  beam 
produced  a  higher  coupling  coefficient  and  Delrin  produced  about  three  times  higher  coupling  coefficient  than  air 
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Figure  14.  Mollier  diagram  for  air  including  singly  ionized  species  derived  from  NASA/GLENN  Chemical 
Equilibrium  Applications  [McBride  and  Gordon  (1996)]  database.  Constant  density  lines  are  indicated  by 
circles  connected  with  dotted  lines.  Molecular  weight  is  indicated  at  intersections  of  isobars  and  isotherms. 
The  1996  database  is  certified  accurate  up  to  20,000  K.  At  20,000  K  air  is  fully  ionized  below  -  /  atm. 
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Figure  15.  Mollier  diagram  for  equilibrium  air  without  ions.  Numbers  adjacent  to  intersections  of  isobars  and 
isotherms  indicate  average  molecular  weight.  Upper  heavy  line  is  a  constant  density  -1.18  kg/m3  line,  three 
isentr  opes  from  u  =  4e7  and  2e7  J/kg  represent  equilibrium  blowdown  for  the  indicated  area  ratios,  1  to  32768. 
The  P  =  1  atm  isobar  is  reached  with  an  area  ratio  -  8  from  the  initial  u  -  4e7  and  ~~  4  from  the  u  =  2e7  J/kg 
state,  respectively.  The  exit  velocities  from  these  two  expansions  are  -  5000  and  3000  m/s,  respectively. 
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Figure  1 6.  Process  representations  of  air  expansion  with  frozen  blowdown ,  as  in  Figure  15.  The  frozen 
isentr  opic  expansions  from  u  -  4e7  and  2e7  J/kg  and  1. 18  kg/m3  to  1  atm  produce  exit  velocities  of 5000  and 
3000  m/s,  only  slightly  less  than  the  equilibrium  expansions  represented  in  Figure  15. 
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Figure  18.  Fractional  energy  conversion  (aft)  as  a  function  of  initial  density  and  area  ratio  for  blowdown  of  air 
heated  to  specific  internal  energy  4e7  J/kg  and  expanded  to  vacuum.  For  initial  density  of  1  kg/m  ,  ct.fi  =  0.50 
for  an  area  ratio  of  4.  afi  values ,  calculated  for  equilibrium  expansion ,  are  only  a  few  %  larger  than  frozen 
ex  vans  ion.  _ _ _ _ 
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Figure  19.  Exit  velocity  from  equilibrium  expansion  of. air  heated  to  u  =  4e7  J/kg  and  initial  densities  ranging  from  le-5 
to  le2  kg/m3  as  a  function  of  area  ratio.  Curves  with  filled  triangles  denote  exit  velocity  that  derives  from  momentum 
thrust  alone  and  are  transformations  of  isentropes  from  the  Mollier  plane  that  originate  from  the  constant  u  =  4e7  J/kg 
line.  Curves  defined  by  squares  show  exit  velocity  for  expansion  into  vacuum  and  include  impulse  from  momentum  and 
pressure  thrust.  Curves  defined  by  circles  show  exit  velocity  for  expansion  to  1  atm  external  pressure.  These  are 
truncated  at  the  point  where  the  impulse  from  pressure  thrust  becomes  negative  due  to  over  expansion  of  the  propellant 
for  the  given  area  ratio.  The  actual  average  exit  velocity  from  the  blowdown  beginning  at  a  given  density  is  the  density 
weighted  average  exit  velocity.  For  example ,  the  heavy  solid  line  represents  the  vacuum  exit  velocity  for  blowdown  from 
an  initial  density  of  1  kg/m3  and  an  area  ratio  of  4.  Initially  the  exit  velocity  is  ~  5800  m/s  and  the  nozzle  exit  density  is  ~ 
0.09  kg/m3.  After  90%  of  the  propellant  has  blown  down ,  the  exit  velocity  is  ~  5500  m/s  and  the  nozzle  exit  density  is  ~ 
0.009  kg/m3.  After  99%  blowdown ,  the  exit  velocity  is  ~  5300  m/s  and  the  nozzle  exit  density  is  ~~  0.0009  kg/m3.  At 
99.9%  blowdown  exit  velocity  is  5100  m/s  and  the  nozzle  exit  density  is  0.0.00009  kg/m3.  Thus  the  average  vacuum  exit 
velocity  is  ~~  5763  m/s ,  is  only  ~  1%  less  than  the  initial  exit  velocity.  In  this  example ,  if  expansion  to  external  pressure 
of  1  atm  pressure  at  area  ratio  of  4  is  considered the  initial  exit  velocity  is  5500  m/s  and  after  90%  blowdown  it  is  ~ 
5000  m/s.  Thus ,  there  is  a  significant  negative  pressure  thrust  involved  when  blowdown  to  1  atm  external  pressure 
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